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We consider the APED model (activation-polymerization-epimerization-depolymerization) for describing the emergence of chiral 
solutions within a non-catalytic framework for chiral polymerization. The minimal APED model for dimerization can lead to 
the spontaneous appearance of chiral oscillations and we describe in detail the nature of these oscillations in the enantiomeric 
excess, and which are the consequence of oscillations of the concentrations of the associated chemical species. 



en 



(D 1 

o. 
x. 

Of 



> 

00 

o 



X 



1 Introduction 

Chemical oscillations are a paradigmatic example of temporal 
self-organization in non-equilibrium systems and have been at 
the center of intense experimental and theoretical research for 
decades 1 " 3 , with theoretical work going back at least to the 
times of Lotka 4 . At the heart of these spontaneous oscilla- 
tions are a positive feedback loop, often an autocatalytic step, 
and a negative control loop, responsible for the saturation of 
the growth process and which operates often on a slower time 
scale. 

The mathematical description of chemical oscillations re- 
lies usually on ordinary differential equations readily derived 
from the chemical rate laws. The number of reacting species 
determines the number of variables and chemical reaction ki- 
netics implies that these equations are nonlinear. The number, 
nature, and stability of the solutions of the system are stud- 
ied through standard methods from dynamical systems theory, 
bifurcation theory, among others^. One simple possibility 
to obtain stable oscillations is the Hopf bifurcation scenario, 
discussed in more detail below. 

Instabilities observed in chemical systems are nowadays 
well-established in theory and experiment, and in tradition 
of the famous Belousov-Zhabotinsky reaction 7 , most exper- 
iments refer to chemical reactions in solution. Nevertheless, 
surface chemical reactions and an increasing number of bio- 
chemical reactions have also been studied extensively 8,9 . 

In this article, we study a model purporting to describe gen- 
eral chiral polymerization processes in a prebiotically relevant 
chemical system, devised for explaining the appearance of ho- 
mochirality 10 . Although there are many different theoretical 
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(and experimental) models tackling the open question of the 
origin of biological homochirality, i.e., the preference of liv- 
ing matter for only one of two chiral states of an otherwise 
chemically identical molecule, at the heart of many of these 
approaches lies a mechanism responsible for spontaneous mir- 
ror symmetry breaking. The initially small chiral fluctuation 
must then be amplified to a state that is useful for biotic evo- 
lution. A precondition for this to happen is that the chirality 
must be preserved and transmitted to the rest of the system. 

In this context, the implications of chiral oscillations for 
chirality transmission are far-reaching: indeed, the memory 
of the sign of the initial primordial chiral fluctuation is erased 
and diluted as it were by whatever subsequent oscillations are 
present in the enantiomeric excess. This phenomenon, if and 
when it occurs in any model purporting to have relevance to 
prebiotic chemistry, adds a further element of uncertainty and 
stochasticity to the determination of the sign of the chirality 
that is finally transmitted to the system (say, within a closed 
reacting system with necessarily damped chiral oscillations). 
Recently in fact, numerical evidence for such damped oscil- 
lations was found in a chiral polymerization model closed to 
matter and energy flow, where the amplitude of the oscilla- 
tion depends on the length of the homochiral chain formed 11 . 
While such models are successful in fitting to real chemical 
data on e.g., relative abundances of chiral polymers 1 -, they 
usually involve dozens or even hundreds of nonlinear coupled 
differential equations, and they do not lend themselves easily 
to a systematic controlled mathematical analysis of the kind 
needed to study rigorously oscillatory phenomena in detail. 

For this essentially technical reason, a detailed study of 
chiral oscillatory phenomena is best carried out in a sim- 
ple model, ideally not too far removed, if possible, in spirit 
from a more complex "real" model. Linear stability analy- 
sis proves that there is no possibility for oscillations in the 
original Frank model, nor within any of its minimal exten- 
sions—. On the other hand, the model proposed by Plasson 
et al.— is amenable to theoretical analysis, and it exhibits dif- 
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ferent stable stationary states, together with instabilities and 
even temporal oscillations. These oscillations refer not only 
to the concentrations of the different species, and therefore 
represent chemical oscillations, but also to oscillations in the 
enantiomeric excess of the system. These chiral oscillations 
are the focus of this article. 

This article is organized as follows: In Sec. [2] we intro- 
duce the minimal APED model and perform a basic bifurca- 
tion analysis. In Sec. [3] we discuss in more detail the basic 
bifurcations observed in the model and display bifurcation di- 
agrams for other parameters. One of the interesting results 
is that chiral oscillations can be expected for parameter val- 
ues that are chemically more realistic, at least within simple 
models based on Frank's scheme, when these are proposed as 
the fundamental candidate reaction network for explaining the 
salient features of the Soai reaction 14 . Finally, we close the 
article with a discussion of the results in Sec. [4] 

2 The APED model 

We consider the so-called APED (activation-polymerization- 
epimerization-depolymerization) model, introduced by Plas- 
son et al.— , used to describe the spontaneous onset of ho- 
mochiral states in a symmetric system of interacting polymer- 
ization products. The reaction scheme for the minimal APED 
dimerization model is represented by the following chemical 
transformations 
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The model contains deactivated monomers L and D, acti- 
vated monomers L* and D*, homodimers LL and DD, and het- 
erodimers LD and DL. The reactions include activation (rate 
a), deactivation (rate b), homochiral polymerization (rate p), 
heterochiral polymerization (rate ap), homochiral hydrolysis 
(rate h), heterochiral hydrolysis (rate j3/i), homochiral epimer- 
ization (rate e), and heterochiral epimerization (rate ye). In 
general, polymerization, hydrolysis, and epimerization pro- 
cesses are stereoselective, i.e., a ^ 1, J3 ^ 1, J ^ 1. Mass 
is conserved and therefore the total concentration in residues 



c = [L] + [D] + [L*] + [D*] + 2([LL] + [DD] + [LD] + [DL]), is 
constant. 

The chemical reactions transcribe into the following set of 
ordinary differential equations: 

[L] = -a[L]-p[L][V] - ap[L][D*] + 2/z[LL]+ 

+ b[L*]+Ph[DL] + plh[LD] (2a) 



[D] = - fl [D] -p[D][D*] - ap[D][L*] +2/z[DD]+ 





+ fr[D*]+j3/i[DL] + /3/![LD] 


(2b) 


[I>] 


= a[L] -p[L][L*] - ap[L*][D] - b[V] 


(2c) 


[D*] 


= a[D] -/7[D][D*] - ap[D*][L] - b[V*] 


(2d) 


[LL] 


= p[L][V] - h[LL] + e[DL] - ye[LL] 


(2e) 


[DD] 


= p[D][D*] - MDD] +e[LD] - ye[DD] 


(2f) 


[DL] 


= ap[L][D*] -e[DL] -j3/i[DL] + ye[LL] 


(2g) 


[LD] 


= (c - [L] - [D] - [L*] - [D*] - 2[LL]- 






-2[DD]-2[DL])/2 


(2h) 



where due to mass conservation, the concentration of one of 
the species, here chosen to be LD, is given by the others. An 
important quantity is the enantiomeric excess in the system. 
Unless otherwise stated, we consider the total enantiomeric 
excess, ee = ([L] + [L*] - [D] - [D*] + 2[LL] - 2[DD]) fc. 

It is straightforward to search for different stable states 
within the system. Following the analysis published in Ref. l0 , 
we present in Fig. Q] a bifurcation diagram obtained with 
MatCont 15 for the a — c space. The different behaviors of 
the system depend on the total concentration c and the rela- 
tive strength of the hetero- vs. homo-dimerization, a. We 
start our analysis at the point c = 1, a = 10 (other parameters 
are a = h = e = p= \ and b = ft = y = 0), where we find a 
stable racemic state (vanishing enantiomeric excess) charac- 
terized by constant and nonvanishing values for all variables 
(but with [D]=[L], [DD]=[LL] etc.). We keep a = 1, b = for 
the remainder of this article. 

If we decrease the total matter in the system (lowering c), 
the system crosses the blue curve and the only stable state is 
the "dead state", characterized by vanishing concentrations for 
all species except for the activated species L* and D*. This 
curve is given by the equation c = 2a/(p + p(X), see Ref.—. If 
we decrease a while c > 1, the system undergoes a pitchfork 
bifurcation at a — 1, at which the racemic state becomes un- 
stable and stable chiral states appear. These chiral states were 
the focus of the article 10 , and therefore we do not discuss them 
in detail here. 

As we increase a for fixed c (which is allowed to have val- 
ues as small as c « 0.2), we observe the onset of oscillations 
as the racemic stationary state becomes unstable through a su- 
percritical Hopf bifurcation (red curve). This means that close 
to the bifurcation, oscillations have small amplitude and fi- 
nite frequency. These oscillations are exemplified in Fig.|2]for 
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Fig. 1 Bifurcation analysis of the APED model in the a — c space for the parameter values a = h = e = p= l and b = p = y = 0. The GH 
point (open circle) is located at GH= (0.19559, 16.23713). Stable oscillations of the variables (and enantiomeric excess) are found in the 
region between the supercritical Hopf curve and the LPC (for: Limit Point bifurcation of Cycles) curve. In (a), a is shown on a linear scale 
from to 60 and the insets qualitatively display the main solutions: open (filled) circle for the unstable (stable) fixed point, i.e., racemic state, 
and solid (dashed) curves for the stable (unstable) limit cycle, i.e., oscillations, (b) Same curves shown on a logarithmic scale for a, for 
comparison with Fig. 3 of Ref. 10 . Note that oscillations cannot be found for any a if c becomes sufficiently large. More information see text. 
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Fig. 2 Stable oscillations of the concentrations and the enantiomeric 
excess ee observed in the APED model for the parameter values 
a = h = e = p = c = l,b = P= y=0, and a = 40. All 
concentrations oscillate, for simplicity we just display [L] and [D], 
and have omitted the initial transient. Note that [L] and [D] oscillate 
out of phase. 



a = 40, c = 1, i.e., already at some distance from the Hopf 
point located at a = a# = 38.24. All concentrations oscillate 
(with a period T = 6.52355), with the chiral species oscillat- 
ing in antiphase, and the enantiomeric excess oscillates with 
the same period. The range where stable chiral oscillations are 
observed does not extend to the hyperbola limiting the dead 
state. As c is decreased for a fixed a, the stable limit cycle un- 
dergoes a saddle-node bifurcation of limit cycles (black curve, 
denoted as LPC, for: Limit Point bifurcation of Cycles). This 
means that the stable limit cycles collides with an unstable 
one and both disappear (and hence no limit cycles are present 
to the left of the black curve). The fixed point describing the 
racemic state remains unstable in that parameter region and 
numerical simulations (not shown here) indicate that the sys- 
tem dynamics is governed by irregular oscillations that do not 
saturate and finally lead to negative values for the concentra- 
tions, hence not representing any realistic chemical state. 

The curve of the supercritical Hopf bifurcation and the 
curve of saddle-node bifurcation of limit cycles meet at a Gen- 
eralized Hopf (GH) point at GH= (0.19559, 16.23713). This 
implies that the supercritical Hopf bifurcation converts into a 
subcritical one. Hence, the red curve between the GH point 
and the point where the Hopf bifurcation meets the hyperbola 
c = 2/(1 + a) represents a subcritical Hopf bifurcation where 
the unstable fixed point (above the curve) becomes stable (be- 
low the curve), representing the stable racemic state, while at 
the same time an unstable limit cycle appears. This unstable 
limit cycle is the one which disappears at the LPC curve. The 
insets of Fig. QJa) show qualitatively the stationary racemic 
state and the limit cycles as one moves around the GH point in 
a — c parameter space. In the parameter regime where stable 



1-1|3 



p 







1 I 1 I 1 


IJ. — 1 — *~ 
— c= 1.0 Hopf 










■■■■ c=1.0LPC 






! 




— c=0.3Hopf 




- 






— c=0.3LPC 


- 


- 


1--/ 


sustained oscillations 


- 




<D / 










Si / 






















1 1 I I I I I I I I I I I I I I I I I I I I 

10 20 30 40 50 60 70 80 90 100 

a 



Fig. 3 Bifurcation analysis in p — a space for the parameter values 
c = 1 and c = 0.3. The area with sustained oscillations is found to 
the right of the Hopf curves and above the LPC curves, respectively. 
The GH and the intersection with c = 2a/ (p + pa) (appearance of 
the dead state) are denoted by open circles and squares, respectively. 
Forc = l: GH= (16.23713,0.19559), intersection at 
(15.36557,0.12221). Fore = 0.3: GH= (16.23713,0.65197), 
intersection at (15.36451,0.40739). Other parameters are 
a = h = e = 1 and b = /3 = y = 0. 

chiral oscillations are observed, another, unstable, limit cycle 
is present. For this parameter set, the minimal a-value for 
oscillations is given by a c = 16.237. 

In Fig. [Tib), we show the same bifurcation diagram in a log- 
linear plot. Contrary to what was displayed in Fig. 3 of Ref. 10 , 
the Hopf curve bends upward for c — > 2, demonstrating that the 
oscillations are not posible if the total amount of residues is 
large. Therefore, the region for c for which stable oscillations 
are found is bound from below and above. 

3 Bifurcation analysis 

The bifurcation analysis carried out in Sec.|2]shows that for the 
parameter values a = h= e = p = \ and b = j3 = y = 0, stable 
oscillations are only to be found for a > a c = 16.237. Here, 
we consider other (obviously not all) regions of the overall 9- 
parameter space (we always keep a = 1, b = 0) to check how 
general is the appearance of chiral oscillations. Since we focus 
on the oscillations, we do not show a complete bifurcation 
diagram. 

In Fig. |3] we investigate for two different total concentra- 
tions c = 1 and c = 0.3, the region of stable oscillations in the 
p — a space. We hence probe at the same time the absolute 
influence of the polymerization rate (p) and the relative influ- 
ence of its stereoselectivity (a). The region of stable oscilla- 
tions is again limited by a Hopf curve and an LPC curve. The 
two curves are qualitatively similar for both values of c, so we 



restrict the description to the case c — 1. If we choose, e.g., 
p = 1, and increase a, the system undergoes a supercritical 
Hopf bifurcation at an = 38.24 as already seen above. Hence, 
to the right of the Hopf curve, the region of stable oscillations 
is found, which is limited from below by the LPC curve. This 
means that if p takes smaller values, smaller a's are necessary 
to induce oscillations, although there is a minimum value for 
p, below which no oscillations can be found, here, p w 0.2, in- 
dependently of a. The LPC and Hopf curves merge at the GH 
point, giving again a minimum value for a = a c , below which 
no oscillations are observed. For smaller c, larger p have to be 
chosen to obtain stable oscillations. 

To measure how the stereoselectivity of polymerization and 
depolymerization influence the appearance of oscillations, we 
perform a bifurcation analysis in j3 — a space (Fig. Ufa)). 
Since c = 1 and other parameters are also as in Fig. [JJ for 
/3 = 0, stable oscillations are found for a > an (at least for 
a < 1500). If < J3 < 0.020445, i.e., for only weak hete- 
rochiral depolymerization, the range of a that allows for stable 
oscillations shrinks, with a lower bound that increases and an 
upper bound that decreases. If /3 > 0.020445, there is no value 
of a that admits stable oscillations. At the same time, from the 
figure we deduce that the range of allowed values of j3 is maxi- 
mal for a = 102.3. Also in this parameter space, we observe a 
Generalized Hopf point at GH= (358.033236,0.011450). To 
the left of that point, the Hopf bifurcation is supercritial, to 
the right it is subcritical. There is a LPC curve emerging at 
the GH point that in this scenario lies above the subcritical 
Hopf curve. This implies that the local bifurcation scenario 
is slightly different to the one displayed in Fig. [JJ The in- 
sets show the main solutions. In particular, we have a stable 
fixed point coinciding with a stable limit cycle between the 
subcritical Hopf curve and the LPC curve. This represents a 
bistability of oscillations with a stationary state. However, this 
scenario is not complete. As we know from Fig.Q] for /3 = 
and a > a.H, there exists also another unstable limit cycle that, 
however, does not undergo any bifurcation in this parameter 
range (and hence is not added to the insets). 

To check the relative impact of the stereoselectivity of poly- 
merization and epimerization, in Fig.Ub), we explore the re- 
gion of stable oscillations in y — a space. Again, the other 
parameters are as in Fig.Q] implying that for y = 0, stable os- 
cillations are found for a > an — 38.24. The functional form 
of the Hopf and LPC curves and the location of GH point are 
qualitatively similar to the scenario in the /3 — a space. There- 
fore, if < y < 0.083136, i.e., for only weak heterochiral 
epimerization, the range of a that allows for stable oscilla- 
tions shrinks, with a lower bound that increases and an upper 
bound that decreases. If y > 0.083136, there is no value of a 
that admits stable oscillations. At the same time, from the fig- 
ure we deduce that the range of allowed values of y is maximal 
for a = 90.05. 
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Fig. 4 (a) Bifurcation analysis in j8 — a space. The area with sustained oscillations is found below the supercritical Hopf curve (to the left of 
GH) and below the LPC curve (to the right of GH). The insets qualitatively display the main solutions: filled (open) circle for the stable 
(unstable) fixed point, i.e., racemic state, and solid (dashed) curves for the stable (unstable) limit cycle, i.e., oscillations. Note that the 
bifurcation scenario is slightly different than in Fig. [7] giving rise to a thin parameter region where bistability among the stable fixed point and 
the stable limit cycle exists. The GH point is located at GH= (358.033236, 0.01 1450), Other parameters are a = c = h = p = e= l and 
6 = y=0. (b) Bifurcation analysis in y — a space. The area with sustained oscillations is found below the supercritical Hopf curve (to the left 
of GH) and below the LPC curve (to the right of GH). The insets qualitatively display the main solutions: filled (open) circle for the stable 
(unstable) fixed point, and solid (dashed) curves for the stable (unstable) limit cycle. The GH point is located at GH= (309.5748,0.043262), 
j8 = 0, and other parameters are as in (a), (c) Bifurcation analysis in /3 — y space. We show three supercritical Hopf curves (a = 50, 100, 200). 
The area with sustained oscillations is found below the corresponding curve. Other parameters are as in (a). 
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Fig. 5 (a) Bifurcation analysis in h — /3 space. We show a supercritical Hopf curve for the parameters a = 50, a = c = p = e = 1, and 
b = y = 0. The area with sustained oscillations is found to the left of the curve, (b) Bifurcation analysis in e — y space. We show a 
supercritical Hopf curve for a = 50, a = c = /? = /i=l, and b = = 0. The area with sustained oscillations is found below the curve. 



To complete the view on relative impact of heterochiral 
polymerization, depolymerization, and epimerization on os- 
cillations, in Fig. Etc), we explore the j3 — y space for fixed a. 
As can be already inferred from the above figures, there is a 
minimum value for a = a# = 38.24, below which no combi- 
nation of (y, j3) yields stable oscillations. However, if a > a#, 
small values of y and j3 are allowed. As we increase a, this 
area first increases and decreases later, once the relative max- 
ima of the /3 — a and y — a curves are passed. 

In an analogous way as our study of p and a, we can com- 
pare the influence of the depolymerization rate h vs. its stere- 
oselectivity j3 (Fig. |5j a)). We choose a = 50 that for j3 = and 
h = 1 (and other parameters as above) corresponds to stable 
oscillations. We find that while for h = 1 the allowed interval 
of j3 is small, we see that if we decrease h, the range of j3 val- 
ues giving rise to oscillations becomes larger. This means that 
the depolymerization process needs not to be very unefficient 
for heterochirals in order to produce oscillations. There exists 
an optimal value h rs 0.2 where the favorable /3 range becomes 
largest. Nevertheless, there is a minimum value h = 0.0931 
below which no stable limit cycles are found. 

Similarly, we study in Fig. |5jb) the influence of the y. In- 
terestingly, we do not observe a minimal threshold for e as for 
h. To the contrary, we find that for e — > 0, the range of allowed 
values of y becomes maximal, and reaches values close to 14. 
However, we should have in mind that if e = 0, the value of y 
becomes irrelevant since ye = 0. Furthermore, since for e = 
the Jacobian matrix becomes singular, a more detailed analy- 
sis should be performed to investigate the bifurcation scenario 
in that limit. We remark that stable oscillations are also possi- 
ble for y=l, i.e., the absence of any stereoselectivity. 

To complete our study, we revisit the j3 — a space, but with 
a parameter set for which 7=1 (Fig. |6). The motivation is 



to find stable oscillations in absence of any stereoselective de- 
and epimerization. We choose c = 1, e = 0.1, and p = 0.3, 
and consider three values for h. We see that the region of 
oscillations increase for smaller h, and in particular, for h = 
0. 1 and /3 = 1, we have a critical value for a = 321 .447 above 
which a stable limit cycle is found. This is just an example to 
demostrate that j3 -C 1 and y -C 1 are not necessary conditions 
for the presence of sustained oscillations. 
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Fig. 6 Bifurcation analysis in /3 — a space. We show three 
supercritical Hopf curves (h = 0.3,0.2,0.1). The area with sustained 
oscillations is found below the corresponding curve. Other 
parameters are a = c = y= 1, e = 0.1, p = 0.3, and b = 0. 
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4 Conclusions 

In this paper, we have studied a simple polymerization model 
(actually, dimerization), aimed originally to explain the emer- 
gence of chiral states, in the context of the onset of chiral os- 
cillations. These oscillations represent temporal changes in 
the enantiomeric excess of the underlying system, i.e., the 
continuous change of a state, where one of the chiral states 
momentarily dominates, into a subsequent state where the op- 
posite handed state is more abundant. 

During decades, models for explaining the emergence and 
amplification homochirality have been focused on the cre- 
ation of stable stationary states, typically associated with chi- 
ral states of maximum enantiomeric excess. The description 
of temporally oscillatory states in these models and the dis- 
cussion of what these chiral oscillations mean, are rare. Chiral 
oscillatory states have been reported by Iwamoto in two dif- 
ferent systems motivated by traditional models for chemical 
oscillations 1617 and shortly after for the APED model 10 stud- 
ied here. Nevertheless, chiral oscillations were not studied in 
detail in the latter paper, and no explanation was given nor rel- 
evance is attributed to this effect. Here, we have studied in 
detail the onset of chiral oscillations, and have demonstrated 
that they are associated with temporal oscillations of the un- 
derlying chemical species, i.e., chemical oscillations. Mathe- 
matically speaking, the onset of oscillations is described by a 
supercritical Hopf bifurcation scenario, implying that oscil- 
lations have small amplitude and finite frequency (opposed 
to other mechanisms where oscillations arise with finite am- 
plitude and vanishing frequency). The Hopf bifurcation is of 
codimension-1, and oscillations are found in large parameter 
regions and hence represent a generic solution of the system. 

While it has been suggested 10 that the reason for observ- 
ing chiral oscillations is that heterochiral polymerization is 
favoured over homochiral polymerization, we test this hypoth- 
esis and clarify in more detail what are the conditions that en- 
able the stabilization of chiral oscillations within the frame- 
work of the rather general APED model. We show continu- 
ation results varying a, j3, 7, p, h, e, and c (keeping a = 1, 
b = 0), demonstrating that oscillations are in fact a typical so- 
lution of the model. In particular, for the a — j3 — 7 space, 
we measure the impact of hetero- vs. homochiral polymer- 
ization, hydrolysis, and epimerization, and find that a> 1, 
P < 1, and 7 < 1 are favorable, in particular a > 38.24 for 
P = 0, and 7=0, where c = 1 is the total mass of the sys- 
tem and other the parameters are also equal to unity. This 
means that oscillations are favored for heterochiral polymer- 
ization and homochiral de-polymerization and epimerization. 
This is a noteworthy result, especially in light of recent anal- 
ysis of Ref. 14 using the Frank scheme to model the Soai re- 
action, which argued that an ideal system for kinetically con- 
trolled spontaneous mirror symmetry breaking should consist 



of a fast exergonic mutual inhibition step leading to the forma- 
tion of the heterodimer LD (i.e., a ^> 1) and which must be 
more exergonic that the possible homodimerization leading to 
LL and DD. The latter condition translates into the constraint 
a 3> ft in terms of the APED model. In all cases observed, the 
oscillations arise from a achiral state. Oscillations are not oc- 
curring for any amout of total mass in the system. We observe 
that there is a lower bound (even before the racemic solution 
reaches the dead state) and also an upper bound. 

We show stable sustained oscillations for a wide range of 
parameter sets. In particular, j3 and 7 can also be chosen to 
unity, i.e., depolymerization and epimerization need not be 
stereoselective in order to produce oscillations, although het- 
erodimers must be polymerized preferentially. Recent exper- 
imental results demonstrate the stereoselectivity of polymer- 
ization reactions and show in detail that the stereoselectivities 
can be modified experimentally (e.g., by pH) in magnitude and 
can be found to be both larger or smaller than unit y 18 ' 19 . 

The presence of oscillations of chirality adds a further el- 
ement of uncertainty to the determination of the sign of the 
chirality that is finally transmitted to the system (say, within a 
closed reacting system with necessarily damped chiral oscil- 
lations). In some parameter regions, even bistability of oscil- 
lations with the racemic state are observed, implying that then 
it is a matter of the initial condition which state is realized. 

Until recently, chiral oscillations have not been reported in 
experiments, and hence the possibility for studying their on- 
set in theoretical models has received sparse attention. How- 
ever, there is now growing interest in oscillatory phenom- 
ena close to the onset of chirality, as chiral oscillations have 
been reported in different experimental systems consisting of 
carboxyclic acids (profens, amino acids, hydroxy acids) 20 - 21 . 
There, several theoretical models are presented and checked 
whether they can account for the experimental observations. 
Among them, the APED model studied here, is the only one 
that was already known to display chiral oscillations, while the 
other ones-2 2 -"— had to be modified or reinterpreted in order to 
account for chiral, rather than just mere chemical oscillations. 
Here, we do not aim to discuss any of those models with re- 
gard to the experiments (the reader is referred to 20 ' 21 directly). 

For the Soai reaction, Buhse presented a simplified model- 
that was later modified to investigate temporal oscillations in 
continuous flow and semibatch conditions 26 . While there are 
stable oscillations of the reactants in a continuous flow sce- 
nario, they are apparently not associated with oscillations in 
the enantiomeric excess. This means that chemical oscilla- 
tions are not related to chiral oscillations, as they are in the 
APED model. On the other hand, in the case of a semibatch 
system, oscillations necessarily fade out and give rise to a chi- 
ral state, but transient oscillations in the enantiomeric excess 
can be observed. 
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